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Abstract 

The recent researches in non equilibrium and far from equilibrium 
systems have been proved to be useful for their applications in differ- 
ent disciplines and many subjects. A general principle to approach 
all these phenomena with a unique method of analysis is required in 
science and engineering: a variational principle would have this fun- 
damental role. Here, the Gouy-Stodola theorem is proposed to be 
this general variational principle, both proving that it satisfies the 
above requirements and relating it to a statistical results on entropy 
production. 



1 Introduction 



The variational methods are very important in mathematics, physics, science 
and engineering because they allow to describe the natural systems by phys- 
ical quantities independently from the frame of reference used p]. Moreover, 
Lagrangian formulation can be used in a variety of physical phenomena and 
a structural analogy between different physical phenomena has been pointed 
out [2] . The most important result of the variational principle consists in ob- 
taining both local and global theories [3]: global theory allows us to obtain 
information directly about the mean values of the physical quantities, while 
the local one about their distribution j3]. 

Last, the notions of entropy and its production and generation are the 
fundamentals of modern thermodynamics and a lot of variational approaches 
has been proposed in thermodynamics [5] . Today, the researches in non equi- 
librium and far from equilibrium systems have been proved to be useful for 
their application in mathematical and theoretical biology, biotechnologies, 
nanotechnologies, ecology, climate changes, energy optimization, thermo- 
economy, phase separation, aging process, system theory and control, pat- 
tern formation, cancer pharmacology, DNA medicine, metabolic engineering, 
chaotic and dynamical systems, all summarized in non linear, dissipative and 
open systems. A general principle to approach all these phenomena with a 
unique method of analysis is required in science and engineering: a varia- 
tional principle would have this fundamental role. But, a useful variational 
principle out of the use of work has been proved not to be obtained J6] and, 
for the open systems, the entropy production for the system and the reservoir 
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has been proved to be obtained if and only if only one heat bath [6] exists and 
its temperature is constant j6]. The different variational principle have never 
been related with these fundamental requirement for a general principle in 
thermodynamics of open system, and often they are related only to closed or 
isolated systems [3]. 

Last, in engineering thermodynamics the open systems are analysed using 
entropy generation, but its statistical definition does not exist and this lack 
does not allow to extend this powerful method to all the above discipline, 
where statistical approach is required. On the other hand, the statistical 
models are not so general to be used in any case with a unique approach; 
consequently, no of them represents a general principle of investigation. 

In this paper, the Gouy-Stodola theorem [7] , [8] , [9] , pH] is proposed to be 
this general variational principle, both proving that it satisfies the above 
requirements and relating it to a statistical results on entropy production. 

2 The open systems 

In this Section the thermodynamic system is defined. To do so, the defini- 
tion of 'system with perfect accessibility', which allows to define both the 
thermodynamic and the dynamical systems, must be considered. 

An open N particles system is considered. Every i— th element of this 
system is located by a position vector x» G -R 3 , it has a velocity Xj G -R 3 , a 
mass mi G R and a momentum p = mji i, with i G [1,N] and p G i? 3 [TTj. 
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The masses rrii must satisfy the condition: 

JV 

\J rrii = m (1) 
i=i 

where m is the total mass which must be a conserved quantity, so it follows: 



where p = dm/dV is the total mass density, with V total volume of the 
system and x B G i? 3 , defined as x B = ^iliP«/ m ; velocity of the centre of 
mass. The mass density must satisfy the following conservation law [I]: 



where p, is the density of the i—th elementary volume with J2i=i = V, 
and S is the source, generated by matter transfer, chemical reactions and 
thermodynamic transformations. This open system can be mathematical 
defined as follows. 

Definition 1. [12] - A dynamical state of N particles can be specified by the 
3N canonical coordinates {q, G R 3 ,i G [1, N]} and their conjugate momenta 
{pi G R 3 ,i G [l,iV]}. The6N— dimensional space spanned by {(pi,Qi) ,i G [l,iV]} 
is called the phase space Q. A point o~i = (pi, q^)^^ ^ in the phase space 
f2 := {(Tj G R 6N : o~i = (p 4 , q^) , i G [1, N]j represents a state of the entire 
N— particle system. 

Definition 2. [11] - A system with perfect accessibility Qpa is a pair (O, 
with II a set whose elements ir are called process generators, together with 



p + pV • ± B = 



(2) 



pi + pi V • ±i = pE 



(3) 
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two functions: 



TT l-> S 



(4) 



TT , TT J I — y TT 7T 



(5) 



where S is the state transformation induced by tt, whose domain T> (tt) and 
range IZ(tt) are non-empty subset ofQ. This assignment of transformation 
to process generators is required to satisfy the following conditions of acces- 
sibility: 

1. Ilcr := {Sa : tt G II, a G T> (n)} = Q , Vcr G Q: the set Ua is called the 
set of the states accessible from a and, consequently, it is the entire 
state space, the phase spase Q; 



2. ifV (n")nll (tt') ^ V (tt'V) = S', 1 (V (tt")) andS^a = S^S^o 
Va G T> (tt'V') 



Definition 3. [11] - A process in Qpa is a pair (tt, a), with a a state and tt 
a process generator such that a is in V (tt). The set of all processes of Sip a 
is denoted by: 



If (tt, a) is a process, then a is called the initial state for (tt, a) and Sa is 
called the final state for (tt, a). 

Definition 4. In an open system, there exists a characteristic time of any 
process, called lifetime t, which represents the time of evolution of the system 
between two stationary states. 



n o Q = {(tt, or) : tt G n, a G V (tt)} 



(6) 
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Any observation of the open irreversible system, in order to evaluate its 
physical quantities related to stationary states, must be done only at the 
initial time of the process and at its lifetime. During this time range the 
system moves through a set of non-equilibrium states by fluctuations on its 
thermodynamic paths. 

Definition 5. W - A thermodynamic system is a system with perfect acces- 
sibility Qpa with two actions W (n, a) G R and H(n,a) G R, called work 
done and heat gained by the system during the process (jr, a), respectively. 

The set of all these stationary states of a system Qpa is called non- 
equilibrium ensemble [TTJ . 

Definition 6. [11] - A thermodynamic path 7 is an oriented piecewise con- 
tinuously differentiable curve in flpA- 

Definition 7. FIT] / - A cycle C is a path whose endpoints coincide. 

Definition 8. WJf - A smooth map S of a compact manifold M. is a map with 
the property that around each point a it can be established a system of coordi- 
nates based on smooth surfaces W* and W™, with s=stable and u=unstable, 
of complementary positive dimension which is: 

1. covariant: dSW^ = W^ a ,i = u,s. This means that the tangent planes 
dSW^, i = u,s to the coordinates surface at a are mapped over the 
corresponding planes at So; 

2. continuous: dSW l a , with i = u,s, depends continuously on a; 

3. transitivity: there is a point in a subsistem of Qpa of zero Liouville 
probability, called attractor, with a dense orbit. 
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A great number of systems satisfies also the hyperbolic condition: the 
lenght of a tangent vector v is amplified by a factor C\ k for k > under 
the action of S~ k if a G Wj* with C > and A < 1. This means that if an 
observer moves with the point a it sees the nearby points moving around it 
as if it was a hyperbolic fixed point. But, in a general approach this 

The experimental observation allows to obtain and measure the macro- 
scopic quantities which are mathematically the consequence of a statistics 
fiE describing the asymptotic behaviour of almost all initial data in perfect 
accessibility phase space Qpa such that, except for a volume zero set of initial 
data cr, it will be: 

1 f 

lim - V up (S J a) = \ [iE (da) up (a) (7) 

for all continuous functions ip on Qpa and for every transformation a i— > 
St (cr). For hyperbolic systems the distribution is the Sinai- Ruelle-Bowen 
distribution, SRB-distribution or SRB-statistics. In particular, here, the 
statistics is referred to a finite time r process, as every real process is, so it 
is considered a SRB-statistics for a finite time system, which exists even if it 
is not so easy to be evaluated. 

The notation iie (do) expresses the possible fractal nature of the support 
of the distribution fiE, and implies that the probability of finding the dynam- 
ical system in the infinitesimal volume da around a may not be proportional 
to da [H]. Consequently, it may not be written as \Le (a) da, but it needs 
to be conventionally expressed as iie (da). The fractal nature of the phase 
space is an issue yet under debate [13] ; but there are a lot of evidence on it 
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in the low dimensional systems |16j . Here this possibility is also considered. 

Definition 9. [11] - The triple (ftp a, J 7 , He) is a measure space, the Kol- 
mogorov probability space, T. 

Definition 10. 07]/ - A dynamical law t& is a group of meausure-preserving 
automorphisms S : Qpa &>pa of the probability space T. 

Definition 11. /IT]/ - A dynamical system T d = (Qpa,F, fiE,T~d) consists of 
a dynamical law r d on the probability space T. 



Irreversibility occurs in all natural processes. In accordance with the second 
law of thermodynamics, irreversibility is the phenomenon which prevents 
from extracting the most possible work from various sources. Consequently, 
it prevents from doing the complete conversion of heat or energy in work; 
indeed, in all the natural processes a part of work, W\ is lost due to ir- 
reversibility. This work can be related to the entropy generation. In this 
Section the entropy generation and its relation to the work lost due to irre- 
versibility is developed for the open systems, introducing the Gouy-Stodola 
Theorem. 

Definition 12. flffl The work lost W\ for irreversibility is defined as: 



3 The Gouy-Stodola Theorem 
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where W\ is the power lost by irreversibility, defined as: 



Wx = W max - W (9) 

with W max maximum work transfer rate (maximum power transferred) , which 
exists only in the ideal limit of reversible operation, and W the effective work 
transfer rate (effective power transferred). 

Definition 13. [TE^ The entropy of the whole system, composed by the open 
system and the environment is defined as: 



$= I \ 5 S- ) =±S e + S g (10) 



T 

where S g is the entropy generation, defined as: 



S g = / dtS g (11) 



with S g entropy generation rate defined as: 

Sg = -7^7 + G ou tS ou t — Gi n Si n — — (12) 

out in i=l 1 

while AS e is defined as the entropy variation that would be obtained ex- 
changing reversibly the same heat and mass fluxes throughout the system 
boundaries, G in the mass flow, the terms out and in means the summation 
over all the inlet and outlet port, s is the specific entropy, S is the entropy, 
Qi,i G [l,iV] is the heat power exchanged with the i—th heat bath and Tj its 
temperature, r is the lifetime of the process which occurs in the open system. 
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Then the term due to irreversibility, the entropy generation S g , measures 
how far the system is from the state that will be attained in a reversible way. 

Theorem 1 (Gouy-Stodola Theorem). In any open system, the work lost 
for irreversibility W\ and the entropy generation S g are related each another 
as: 

W x = T a S g (13) 
where T a is the ambient temperature. 

Proof. Considering the First and Second Law of Thermodynamics for the 
open systems, the maximum power transferred is: 

W max = Y,G l nU+ V -+gz+T a s\ -Y, G outU+ V -+gz+T a s\ -j(E-T a S) 

in ^ " ' in out ^ ~ ' out 

(14) 

while the effective power transferred results: 

W = Y,G m z+T a s) -J2 °out ( h +t+9 z+T a s) -J,(E-T a S)-T a S g 

in \ J in out V ~ / out 

(15) 

where h is the specific enthalpy, v the velocity, g the gravity constant, z the 
height and E is the instantaneous system energy integrated over the control 
volume. 

Considering the definition of power lost, W\, it follows that: 

W x = T a S g (16) 
form which, integrating over the range of lifetime of the process, the Gouy- 
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Stodola theorem is proven: 



/ 

JO 



T 



J 

Jo 



T 



dtW; 



T, 



a 




(17) 



□ 



The Gouy-Stodola result on the entropy generation is expressed in a 
global way, without any statistical approach and expression. In order to 
extend the use the Gouy-Stodola theorem to any approach and context, a 
statistical expression of entropy generation is required. To do so the following 
definition can be introduced: 

Definition 14. J7Pj / The entropy production Yi wol i is defined as: 



with N number of heat baths, whose temperature is Ti, i £ [1, N] in contact 



k,B Boltzmann constant, £(er) phase space contraction and SRB- statistics. 

Theorem 2. In a stationary state, the entropy generation and the entropy 
production are related one another by the relation: 




N 




with the system, Qi, i £ [1, N] heat power exchanged with each i—th heat bath, 




(19) 



Proof. If the system is in a stationary state (dS/dt = 0) and if: 



1. the system is closed: Y^out G 0U tS ut 



and GinSin, or 
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2. the system is open, but Y.out G 0U tS ut - J2in G inSin = 



and considering the relations (12) and (18), then it follows: 



S g — —k B S pr0(j (20) 

from which, integrating the statement is obtained. □ 
Theorem 3. The thermodynamic Lagrangian is: 

C = W X (21) 

Proof. Considering the entropy density per unit time and temperature ps, 
the Lagrangian density per unit time and temperature pc, the power density 
per unit temperature p^, and the dissipation function 0, the following relation 
has been proven [2D] : 

pc = Ps ~ pw - (22) 
and considering that ps — p n = 20, it follows |4J that: 

Pc = (23) 

consequently [TT] . 

£ = fdt [ dT I VpcdV = fdt [ dT f V Pc <pdV = W x (24) 
Jt Jt J Jt Jt J 

□ 

Theorem 4. At the stationary state, the work A lost for irreversibility is an 
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extremum. 



Proof. Considering the definition of action A it follows that: 



A= f dtC= f 
Jo Jo 



dtW; 



A 



(25) 



and considering the least action principle it follows: 



5 A < 5W X < 



(26) 



which allow to state that: 

1. W\ is minimum if the work lost is evaluated inside the system 

2. W\ is maximum if the work lost is evaluated outside the system, inside 
the environment 



Of course, this extremum is extended also the entropy generation, S g , 
using the Gouy-Stodola Theorem W\ = T a S g . 

4 Conclusions 

A general principle to approach the stability of the stationary states of the 
open systems is required in science and engineering because it would repre- 
sent a new approach to the analysis of these systems, with the result of im- 
proving their applications in mathematical and theoretical biology, biotech- 
nologies, nanotechnologies, ecology, climate changes, energy optimization, 



in accordance with the thermodynamic sign convention. 



□ 
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thermo-economy, phase separation, aging process, system theory and con- 
trol, pattern formation, cancer pharmacology, DNA medicine, metabolic en- 
gineering, chaotic and dynamical systems, etc.. 

Here, the Gouy-Stodola theorem has been proved to be the searched 
variational principle, which satisfies the two fundamental request: 

1. to be a work principle, because a useful variational principle out of the 
use of work has been proved not to be obtained [6] 

2. to use only one temperature which remains constant (because the en- 
vironmental temperature is always considered constant in the usual 
applications), because for the open systems, the entropy production 
for the system and the reservoir has been proved to be obtained if and 
only if only one heat bath [B] exists and its temperature is constant [B] . 
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